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Abstract 

We consider the baryon number generation by charge transport mechanism in 
the electroweak phase transition taking properly into account thermal fluxes through 
the wall separating true and false vacuum in the spatial space. We show that the 
diffusion from the true vacuum to the false one has a large diminishing effect on the 
baryon number unless the wall velocity is near to, but less than, the speed of sound 
in the medium and the ratio between the collision rate and wall thickness is about 
0.3. The maximum net baryon density generated is Pb/s ~ 0.2 x 10~ 10 , where s is 
the entropy density of the Universe. If the wall proceeds as a detonation, no baryon 
number is produced. 



It was pointed out by Sakharov ]I] that the cosmological asymmetry between 
baryons and antibaryons can follow from the particle interactions in the early Uni- 
verse if the following conditions are simultaneously satisfied: the baryon number B 
is not conserved, C and CP are not exact symmetries and the Universe is out of 
thermal equilibrium. The electroweak theory has all these ingredients, if the weak 
phase transition is first order ||. The baryon number non- conservation follows from 
the anomalous violation of the B + L number (L is the lepton number) due to 
non-perturbative sphaleron effects . In a first order phase transition there will be 
local departure from thermal equilibrium in the vicinity of the walls between the 
symmetric and broken phases. In the standard electroweak model there also is built 
in CP violation manifested as a complex phase in the Kobayashi-Maskawa quark 
mixing matrix. The resulting CP violating effect is, however, at least two orders of 
magnitude too small to account for the observed baryon to entropy ratio ~ 10~ 10 . 
One thus have to go beyond the minimal standard model. 

One possible mechanism, so-called charge transport mechanism, that may pro- 
duce a large enough baryon asymmetry was suggested by Cohen, Kaplan and Nelson 
PJ. The idea is that a net charge, e.g. lepton number, is reflected from the domain 
wall between broken and unbroken phase, and anomalous SU(2) effects in the un- 
broken phase transmute this into a net baryon number. Then the baryons produced 
pass thermally to the broken phase where they are stable. In order to produce a net 
charge, the reflected particles have to undergo CP violating interactions in the wall. 
In J!] these are the Yukawa interactions of neutrinos which may involve unremovable 
phases. Such situation may be realized, for example, in the context of the singlet 
Majoron model and the left-right symmetric model ||. 

In this paper we shall reconsider some aspects of the charge transport mechanism. 
We derive coupled Boltzmann equations for the evolution of the baryon and lepton 
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numbers in the broken and unbroken phases. In particular, we shall take into account 
also the diffusion of particles from one phase to another, and we actually find out 
that the diffusion has in many cases quite large diminishing effect on the creation 
of the baryon number. 

We are considering a system of an expanding bubble of the true vacuum, i.e. of 
the broken phase. The expansion starts with a critical bubble which is a spherical 
configuration of the order parameter . We assume that the bubble remains spher- 
ical during the expansion, its radius being R(t) = R + v w t, where v w is the velocity 
of the moving wall and Ro is the radius of the critical bubble. The time needed for 
completing the phase transition depends on the size of the critical bubble. In the 
following we will assume a weak phase transition in which case the critical bubble is 
large and the growth time of the bubble is of the same order than the bubble size. 
The growth time is given by M 



0.073 Mpi 

^growth ~ 



ln M « 



7^4 

c 



-3/2 

(1) 



where g* is the effective number of degrees of freedom (here g* ~ 110), Mp\ is the 
Planck mass and T c is the critical temperature. Using the result of || one finds out 
that the bubble radius reaches the value 

R 1 = (8ir)V 3 v w t growth ~ 3.2 x 10V/GeV (2) 
upon completion of the phase transition. The critical bubble radius is thus 



R = Ri- v w t growt h - 2.1 x 10V,/GeV. (3) 

Strictly speaking, the results above are only for the Standard Model, but they 
are likely to give the correct order of magnitudes for a wide range of models not 
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too different from the Standard Model, like the Majoron model. Moreover, our final 
results are quite insensitive to the precise values of Rq and R\. It turns out that 
the changes on R$ and R\ between as large values as 10 4 f w /T and 10 u v w /T do not 
effect the value of the generated baryon number but by a factor 2, at most. 

In the case of a weak first order phase transition nucleated bubbles grow as 
deflagrations. The advancing phase boundary is preceded by a supersonic shock 
front which hits the fluid in the false vacuum and makes it to move to the same 
direction as the boundary. The shock front discontinuity is formed within a time 
period of to = O(10) 1/T C which is very short compared with t growt h @- The velocity 



of the shock front can be estimated as 10 



V shock 



3T , 4 + T /^ 1 (4) 



since the difference between the temperature of the symmetric phase Tf and the 
temperature in the region between phase and shock fronts T q is estimated to be 
smaller than one percent §J. Hence the shock front moves essentially with the 
speed of sound. In front of the shock wall the particles are at rest in average; there 
is no net flux in this region. Inside the shock wave particles are moving with a slow 
velocity outwards causing a uniform particle flux. Behind the wall in the broken 
phase particles are again at rest in average. 

Let us now derive the evolution equations for the lepton numbers L\ and L 2 , and 
for the baryon numbers B\ and B2, where the subscript 1 refers to the symmetric 
region and the subscript 2 to the broken one. We start by writing the equation of 
continuity for the lepton and baryon number densities, Zj and 6j, and the radial fluxes, 
Aj and fa, in the symmetric and broken regions. Assuming spherical symmetry the 
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equations read 

^1 + —^) = -^-[(l-x)b 1 + (l+x)h), 

9 ^ + ^(r%) = 0. (5) 

The source terms in the symmetric region are due to anomalous B + L-violating 
processes, characterized by a time scale tb — lGeV -1 ||. The parameter x is 
defined as the equilibrium value of the ratio (B + L)/(B — L). In the broken region 
the source terms can be neglected, since there the B+L violation effects are expected 
to be small []. 

We now integrate the first two equations of (|5]) over the volume between the phase 
boundary with a radius R(t) and the shock front with a radius R s (t) = R s o + v s t. 
Knowing the exact value of the initial shock front radius R s q is not that important 
as the shock wave is formed, as mentioned, in a very short time period in comparison 
with the whole span of the time evolution. 

The ensuing equation for the lepton number L\ in the symmetric phase reads (L 
stands for time derivative) 



Li = 4irR%(R s ,t)v s - 4vr J R 2 (t)/ 1 ( J R, t)v w + 4vr J R 2 (t)Ai(i?, t) 

- J-[(l- x )B 1 + (l + x)L 1 ]. (6) 
Zt b 

An analogous equation can be derived for the baryon number B\. In deriving Eq. 

(H) it was assumed that the lepton number flux at the shock front can be neglected. 

1 We assume that the sphaleron energy is large enough to prevent the baryon number 
dilution in the broken phase. 
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From the latter two equations of (|5|) one obtains similarly the time evolution 
equation for the lepton number in the broken phase: 

L 2 = 4ttR%{R, t)v w - AtcR 2 X 2 {R, t), (7) 

and an analogous equation for the B 2 . In deriving the equations for L 2 and B 2 it 
was assumed that the flux of the lepton number is negligible at r = Rq for all times. 
Strictly speaking this is not true but serves as a good approximation. 

In order to solve Eqs. (|6|) and (0) and the corresponding equations for the baryon 
number, one has to know the number densities and the fluxes of leptons and baryons 
at the boundary of the two phases. There are different effects one should take into 
account. The CP violating scattering of neutrinos in the advancing wall produce a 
net lepton number flux on the both sides of the wall. This effect was inspected in 
detail in 0] . They found numerically that the thermally averaged net lepton number 
flux Jl, in units of entropy density, reflected from the wall to the symmetric phase 
is 



— fi, = 2.0 x 10 — 4.2 x 10 (8) 

at T = 200 GeV, depending on the wall velocity and thickness. 

There is also a flux of leptons and baryons due to thermal transport of particles 
from one phase to another. The driving force of the thermal flux is the particle 
density differences over the phase boundary. The thermal flux of the lepton number 
in the symmetric phase, to be denoted by Ai, can be presented in the form Ai = 
— DVIj. Here V/i is the radial lepton number gradient and D is the diffusion 
coefficient. A simple estimate for the diffusion coefficient is D = (v) 2 r/3, where 
r is the average time between particle collisions [Q. Because all particles in the 



unbroken phase are massless and they are still relativistic in the broken phase (except 
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the heavy right-handed neutrinos) the average particle velocity (v) is about unity. 
The collision time r can be calculated the model given. In the Standard Model one 
can estimate r ~ O.lGeV -1 . 

If the wall between the broken and unbroken phase were at rest we could ap- 
proximate thermal flux at the boundary as D(li(R(t),t) — l2(R(t),t))/5, where 5 is 
the thickness of the wall. The wall is, however, moving with velocity v w . Taking 
that into account the formula modifies to 

A i = 5 ft W) ,t)(l + v w )-l 2 (R(t) ,t){l-v w )\. (9) 
2d 

The factors (1 + v w )/2 and (1 — v w )/2 give the fractions of light particles which 
are caught by the wall in the symmetric region or catch the wall in the broken 
region, respectively. Those low-energy leptons, which reflect back to the symmet- 
ric region because of becoming too heavy to penetrate the wall, produce a flux 
2r u Dli(R(t), t)(l + v w )/{25). Here the quantity r v gives the fraction of the heavy 
neutrinos that are reflected, and the overall factor of two is there due to the fact 
that in the reflection the lepton number changes by two units. In our numerical 
calculations we will use r„ = 1/4 corresponding to the situation in the Majoron 
model 0. Note that there is no contribution of heavy leptons coming from the 
broken region to the symmetric region, because they are, in average, much slower 
than the wall. The thermal fluxes for the baryon number obtained similarly. We 
have approximated the collision time r to be the same for the both cases. 

About the net lepton and baryon number densities we assume a uniform dis- 
tribution inside the expanded bubble and a linearly decreasing densities inside the 
shock front, i.e. between the radii R(t) and Ro, so that h(R s ) = bi(R s ) = 0. This 
yields 
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L x = -h(R,t)(R s - R)(R 2 S + 2RR S + 3R 2 ), 

O 

4-7T 

L 2 = -jl 2 (R,t)(R 3 -R 3 ). (10) 

Using these approximations the evolution of the lepton and baryon numbers is 
described by the following set of equations: 



= l--^—[(l + v w )(l + r u )a 1 (u)L 1 -(l-v w )(l-r u )a 2 (u)L 2 \-a 1 L 1 



3u 2 du 6Sv 

Rq 



2v w t b 

1 dB 1 



[(l-x)B 1 + (l + x)L 1 ], 



3u 2 du 65v w 

Rq 



f 

— [(1 + v w )ai(u)B 1 - (1 - v w )a 2 (u)B 2 ] - ol\B x 
1 - x)B 1 + (1 + x)Li] , 



2v w t b 

^2~^ = -i + -^-[(^ + v w )(a 1 (u)L 1 -(l-v w )a 2 (u)L 2 \(l-r l ,) + a 2 L 2 , 
3^^7 = ^[( 1 +^)( a i(«) B i-( 1 - | '»h(«) 5 2] + a252. (11) 

The lepton numbers Li and baryon numbers B iy i — 1,2, are presented in units of 
f L and they are scaled with a reference volume V = 47r/3i?Q. We have defined a 
dimensionless variable u = R(t)/Ro and the functions 



ai(u) = A[(u - u)(u 2 + 2uu + 3u 2 )} \ 

a 2 (u) = l/(u 3 -l), (12) 

where u = R s /Ro = R s o/ Ro + v s /v w (u — 1). The baryon density in the broken phase, 
the quantity we are interested in, is then given by 



b 2 = h^B 2 (u = R 1 /R ). (13) 
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Before moving to the numerical solutions of our final equations (|TID , we will 
make some general observations. It is quite easy to deduce from the equations that 
there are solutions with no baryon number generation. This depends on the values 
of the unknown parameters like the wall thickness 5 and the wall velocity v w , as well 
as the diffusion time scale r. It turns out that the decisive quantity is 

Q = <k {1 - vM1 - r ' y (14) 

This can be seen as follows. First make the ansatz that the solution is "almost 
trivial", i.e. B±, B 2 and L\ are zero. Then L 2 of the following form solves (|TT|): 

L 2 = ±-(u 3 -l) + C(u 3 -l) 1 - Q . (15) 

When Q > 1 the initial condition L 2 (u = 1) = forces the integration constant C 
equal to zero. By substituting the solution (|i~5D to eqs. ( p]) one then finds that, 
as consistency requires, the derivatives of L±, B\ and B2 vanish. In order to have 
baryon number generated one has therefore to require Q < 1 . Physically this means 
that the wall have to be thick enough or it must advance into the unbroken phase 
rapidly enough to avoid dilution of the net baryon number due to particle diffusion 
through the wall. Note that this is also the asymptotic solution for the Boltzmann 
equations when R\/Rq ^> 10 9 . Physically that is, however, not a favoured case. 

Let us now present our numerical results. The theoretical estimates for the values 
of the wall velocity v w and the wall thickness 5 vary with the assumptions made of 
the strength of the phase transition and the nature of the bubble expansion. We 
keep these quantities as free parameters. As we have found, the wall thickness in fact 
appears in our expressions always along with the free time r, the results depending 
on their ratio r/8 only. 

In Fig. 1 we present the values of the generated baryon density b 2 as a function 
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of v w and r/8. The maximum value of 62 achieved inside the 62 = 0.2 contour is 
about 0.23, and it corresponds to the wall velocity close to the velocity of sound and 
the wall thickness close to, but slightly larger than, the neutrino free path. We have 
plotted also the Q = 1 contour, to right of which no baryon number is generated, 
because there the wall moves so slowly or is so thin that the particle diffusion through 
the wall dilutes the baryon number. On the other hand, the decreasing of the net 
baryon number with decreasing r/5 follows from the fact that in the case of a thick 
wall less baryon number moves into the broken phase from the unbroken phase where 
it is created. This decrease is of course quite independent on the wall velocity as 
the baryon number flux goes towards the advancing wall. 

In Fig. 1 we have integrated the evolution equations from the initial bubble 
radius Rq to the final radius Ri, whose values were given in eqs. (0) and (|j). As 
mentioned, this corresponds to a weak phase transition where the initial bubbles are 
quite large and the expansion time small. We have also investigated the sensitivity 
of our result on this assumption. This may be relevant since it is also possible that 
the initial bubbles are very small and that the bubble radius increases several orders 
of magnitude upon the completion of the phase transition 0. In Fig. 2 we present 
the generated baryon number as a function of the ratio Ri/Rq with R\ fixed to the 
value (0). One can conclude from this figure that the net baryon number does not 
depend significantly on the initial size of the bubble in the physical region of the 
ratio Ri/Rq. 

We have presented our results in terms of the reflected lepton number flux f^. 
This quantity has been estimated in J| for some values of the parameters. For 
8 ~ 0.1 GeV -1 , v w = 0.87 and the critical temperature T c = 200GeV they give 
fi = —1.2 x 10~ 10 s, where s is the present entropy density in the Universe. Applying 
this value to our results would yield the maximum net baryon number density Pb/s — 
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0.24 x 10 -10 . This is quite close but slightly below the allowed range Pb/s = (0.4 — 
1.0) x 10~ 10 obtained from the standard nucleosynthesis [[0|. The maximum value 
of Pb/s corresponds the values r/8 ~ 0.3 and v w = 0.55. It is relatively sensitive to 
the actual values of wall velocity and ratio r/8. The present theoretical estimates 
for these are 5 ~ (10 — 40)/T c and v w lies between 0.2 and 0.8 ||. Note, however, 
that if the bubble wall proceeds as a detonation, i.e. v w > v s , no baryon number 
is created. In that case there is no region preceeding the bubble wall where lepton 
number conversion to the baryon number could take place. 

Our result seems to give some indication that the baryon number production 
by charge transport mechanism may not be powerful enough to create the observed 
baryon number. Because there is some uncertainty in the actual value of fi, the 
baryogenesis via this scenario is not completely ruled out. The value of fi has been 
calculated only for two wall velocities and for a few wall thicknesses using critical 
temperature T c = 200 GeV. Also the baryon number production might be enhanced 
if there existed relatively large initial baryon and/or lepton number density created 
before electroweak phase transition. Such could be created in the Majoron model by 
decays of of right-handed neutrinos jnj or, in e.g. left-right symmetric model, by 
some mechanism during the symmetry breaking of the right-handed sector. These 
issues are currently under study [O . 
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FIGURE CAPTION 
Figure 1. The net generated baryon density in the (log(r/<5), v w ) -plane in units of 
—jh- With the parameter values corresponding to the area below the Q = 1 curve 
no baryon number is generated. v s is the sound velocity. 

Figure 2. Dependence of the generated baryon density on the ratio Ri/Rq with 
r/6 — 0.2, v w = 0.5 and R\ fixed to the value (0). 
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